Random Schrodinger operators on long boxes, 
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Abstract 

It is conjectured that the eigenvalues of random Schrodinger operators at the lo- 
calization transition in dimensions d > 2 behave like the eigenvalues of the Gaussian 
Orthogonal Ensemble (GOE). We study a the eigenvalues of long boxes in dimension 
d = 2 for low disorder. We deduce a stochastic differential equation representation for 
the limiting process. We show that in dimension d = 2 there are sequences of boxes so 
that the eigenvalues in low disorder converge to Sinei, the limiting eigenvalue process 
of the GOE. 

1 Introduction 



When 



Wignerl ( 119571 ) introduced random matrices to model large atomic nuclei, his main 



goal was to find a simple mathematical model that shows repulsion between eigenvalues as 
observed in the data. The Gaussian orthogonal ensemble (GOE) has since been a remarkable 
success in the physics literature: several important complex systems are predicted to have 
eig envalue repulsion akin to the GOE, most notably the Laplacians of many planar domains, 



sec 



Bohigas et all fll984f ). 



Given the all the non-rigorous theories that predict GOE behavior, as well as a lot of 
numerical evidence, it is surprising that there are hardly any mathematically rigorous results 
in this direction. Most objects with rigorously known bulk GOE behavior are themselves 
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rando m ma trices constr u cted si milarly to the Gaussian orthogonal ensemble, see lTao and Vu 
( 2009 ) and lErdos et al. ( 2009a ) for recent breakthrough results in this direction. 

The goal of this paper is to show that random Schrodinger operators on long boxes in Z 2 
behave like the Gaussian orthogonal ensemble. More precisely, we will show 

Theorem 1. There exists a sequence of weighted boxes on Z 2 with diameter converging to 
oo so that the rescaled eigenvalue process of the adjacency matrix plus suitable diagonal noise 
converges to the Sine\ point process. 

Recall that the Sinei point process is the li mit of the eige nvalues at the bulk of the spec- 
trum of the Gaussian orthogonal ensemble, see Mehtal ( 2004 ). The process-level convergence 
described here implies convergence of eigenvalue gap sizes and all similar local st atistics. 



The concept of localization for random Schrodinger operators was introduced by 



Anderson 



( 119581 ). As more and more diagonal noise is added to the Laplacian of a large box, the eigen- 
vectors change from being sp read out over the entire box to l ocalized on sm aller regions. It 
is believed (see, for example, Altshuler and Shklovskil ( 1986 ). 



Efetov 



(119971 )) that near this 



transition the eigenvalue process is GOE-like, and, as the noise increases, it is approximately 
Poisson. This transition from GOE to Poisson can be analyzed in the long-box case. We 
have 

Theorem 2 (Eigenvalue process transition). Consider boxes in Z 2 with fixed base Z m . Con- 
sider the process of eigenvalues of 

n(l-Z m xZ n + ^=V -IX*), 
2 Jn 



where V are diagonal matrices with independent, mean zero, variance 1 entries with bounded 
third moment. Then 

• For almost all A* G [0,1], for appropriate subsequences, this process has a limit A a 
which depends on m, A*. 

• As o — y 0, the the process A a converges to the (randomly shifted) eigenvalue process of 
a random matrix with independent Gaussian entries. 

• For every a, the process A a is the zero process of the determinant of a m+1- dimensional 
matrix-valued analytic function described by a stochastic differential equation. 

Remark 3. To illustrate the GOE-to-Poisson transition, it is possible to show that 

• as a oo the process A CT converges to the Poisson point process. 
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We plan to do this in a subsequent article. 

The proofs of both theorems are based on the analysis of transfer matrices. The more 
precise statements can be found in Proposition HI Corollary [6], Proposition HI and Corollary 

The noise term in the t ransfer matrix evolution in t his sc aling regime for block Jacobi ma- 



trices has been studied in 



Romer and Schulz-Baldes 



(hoiol ) 



using a language different from 



SDEs. The first arxiv vers i on of the present paper was followed by the preprint of the paper 

who, in independent work, also study SDE limits of trans- 



Bachmann and De Roeck 



fer matrices. Their starting point the so-called DMPK theory in the physics literature, which 
is essentially the study of diffusive limits of quasi-one- dimensional random Schrodinger oper 
ators from a slightly different point of view. We refer the reader to 



Bachmann and De Roeck 



( 120101 ) for a discussion of this theory. One of the novelties of our approach is that it allows 
for studying the dependence on the eigenvalue A, which in turn allows us to deduce the 
scaling limit of the spectrum, the main focus here. 

Noise explosion. The central reason for the appearance of random matrices is noise ex- 
plosion. For a simple case of this, consider 

\nt\ 

where the Xk are independent random variables, with, say, standard normal distribution. 
If 77 = (mod 2tt), then the sequence of functions B v>n converges in law to real Brownian 
motion. Otherwise, B v>n converges to complex Brownian motion. Thus a one-dimensional 
noise process Xk gives rise to a two real-dimensional Brownian motion in the presence of 
oscillatory terms. Moreover, if 77, v are linearly independent over the integers, then B vn , B vn 
converge jointly to two independent Brownian motions. 

Now we change the setting so that the Xk are mxm diagonal matrices with independent 
standard normal entries. Let U be a unitary matrix with eigenangles that are linearly 
independent over the integers, and also assume that the absolute squares of the eigenvectors 
of U are not orthogonal to each other. Now consider the sum 

\nt\ 

A simple computation shows that B n (t) converges to a Hermitian matrix- valued Brownian 
motion process which is m 2 -real dimensional. In such oscillatory sums, noise that was origi- 
nally restricted to the diagonal explodes into all possible directions, and changes dimension 
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from m to m 2 . This phenomenon, which we call noise explosion, plays a central role in 
the proofs below. 



The method of analysis in th is paper is an e x tensio n of the 1-dimensional case in- 
tra duced_jn_j3 L _diffe^^ (120091 ) and further refined and simpli- 
fied in iKritchevski. Valko. and Viragj ( 120111 ). The latter paper studies the most natural 



1-dimensional random Schrodinger operator at the localization transition. 

The current paper does not give any bounds on the aspect ratio of the long boxes m/n. 
This would require rates of convergence estimates for diffusion approximation. 

Problem 1. Show that the result in this paper holds for n = 0(m q ) for some q > 1. 
What is the smallest possible value of q for which the results hold? Is it q = 1? 

Problem 2. Extend the results to long boxes of higher dimension. The eigenvalue 
structure of the base here is more complicated, see Remark [T2j 

Problem 3. Eliminate the r — > condit i on in Proposition [9j This would require more 



general versions of the results in 



Erdos et al. 



2 Description of the model and notation 

Let Z n denote the graph of a path of length n with vertices 1,2, ... ,n. We will use the 
same variables for a (weighted) graph and its adjacency matrix, so Z n will also denote the 
n x n matrix with entries (% n )i,j — l{|i-j|=i}- Let x denote the cartesian product of weighted 
graphs (or matrices), namely 

[A x B)^^^^ = li> = j'A it j + l i= jBi/j> 

For a positive real r and weighted graph (or matrix) G we will denote by rG the graph where 
the weights are all multiplied by r (which is the same as multiplying the matrix G by r). 
Boxes. The boxes we consider are [r(Z m )] x Z n for some positive real r, where n is typically 
much larger than all the ??v This is the adjacency matrix of a weighted graph on the box 
Z m x 7L n where the edges in the first direction have weight r, and the edges in the second 
direction have weight 1. 

The adjacency matrix corresponds to Dirichlet boundary conditions. For a probabilistic 
interpretation, the matrix 

rZ m x Z„ — 2(r + l)I m xn 

is the transition probability matrix of the continuous time random walk on the box killed 
when it leaves; the jump rates are r and 1 in the two different kinds of directions. 
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General slabs. Some of our results will apply to rG x Z n where G is a general symmetric 
matrix. 

Chaoticity. For the noise explosion phenomenon to work, we need to consider conjugation 
by special angles. Consider points (angles) X = {x±, . . . ,x m } on the unit circle M./2n. We 
will be interested in the set 

^ je} U {x h + x i2 + x i3 - x h } 

where the indices run through 1,. . . , m and the operations are meant mod 2n. The chaotic- 
ity of X is defined as the distance of the set A from 

cha(^) = dist(0, A) (1) 

We call the Xi chaotic if they have nonzero chaoticity. 

In particular, if the X are chaotic then then the x^'s are distinct. Moreover, if Xi s are 
linearly independent over the integers then t hey are chaoti c . A condition very similar to 
chaoticity appears as the main assumption in 



Schulz-Baldesl (12004h . 



The noise explosion described in the introduction will rely on the following quantities. 
Let G = ODO^ 1 be the diagonalization of G, so that O is orthogonal and D is diagonal. 
Let Oi = Oi t . denote the ith row of O. Let \Oi\ 2 denote the vector formed by the squares of 
the entries of Oi. 

When G is the adjacency matrix of the discrete interval, then the eigenvectors are of the 
form 



O jk = y/2/(m +l)sm(ir jk/(m + 1)), 

and it is easy to check that 

3/2 for i—j, 

(m + l)(|0. t | 2 ,|0,| 2 )= (2) 

I 1 for % 7^ j . 

3 The regularized transfer matrix evolution 

For this section, consider the operator GxZ„ where G is a symmetric matrix of dimension 
m x m. 

Suppose that we want to solve the equation Mu = \u where M = G x Z n + V with a 
diagonal matrix V. We will do this by solving the system of linear equations recursively, 
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slice by slice (where a slice is a copy of G, i.e. G x {i}). Let u be a function from the vertices 
of G x Z n to M, and let = w(j, i) so that is a vector indexed by G. Denote by Vi the 
diagonal matrix which is the restriction of V to the indices G x {z}. Then the eigenvalue 
equation for entries with the same first coordinate 2 < i < n — 1 reads 



\Ui = (G + Vi)Ui + m-x + u i+1 



(3) 



and if we set no = u n+ \ = 0, then ([3]) holds for i — 1 and n as well. It is not hard to check 
that this system of n equations is then equivalent to Mu = Xu. Note that ([3]) can be written 
in terms of block transfer matrices of dimension 2m as follows: 



Uk+l 









\I-G-Vk 
I 



Uk 

u k -i 



Denoting the above matrix by T^, we can now characterize the eigenvalues of M as 



{A : 3ui, u n such that T n ■ ■ -T 2 Ti 



Ml 















} = {A : det((T n ---T 2 T 1 ) {1 , 1) )}. (4) 



Here the subscript (1,1) refers to the top left m x m submatrix, and both representations 
are equivalent to this submatrix having some eigenvector u\ with eigenvalue 0. So we may 
write 



spec(M) = zeros A [det ((T n • • • T 2 T 1 ) {1)1) )] 



(5) 



where spec(-) will refer to the eigenvalue counting measure and zerosA refers to the zero 
counting measure of a real-analytic function of A. It is not hard to check that the equality 
also holds in the presence of multiple zeros. 

We now introduce the unperturbed version of the transfer matrix 



T 



XJ-G -I 
I 



where A* is some fixed reference point. Let G = ODO T be the diagonalization of G with 
O orthogonal and D = diag(di, . . . , d m ) diagonal. We will use the shorthanded notation 
X Y = Y~ l XY which means G = D° T . 

We may change basis to study products of 



rpO®I 2 
± k 



XI 



D 
I 



0®I 2 



O 

o 
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and note that 

r XJ-D -I 
I 

This matrix can be completely diagonalized as it is just the permutation of a block diagonal 
matrix with 2x2 blocks. (We ha ve learned that such a change of basis has been considered 
for a slightly different problem by ISchulz-BaldesI (120041 ) .) Let 



Zj = (A* - dj)/2 + i\Jl - ((A* - dj)/2y 



and Z = diag(zi, . . . , z m ), We assume that |A* — dj\ G (0, 2) for all j which means that Zj is 
a unit length complex number in the upper half plane. 
Then we have 



where 



and 



S 



diag \ zj 



3 z j\ 







rp(0®I 2 )Q __ 

'1/2 

diag \zj 
Q- 1 = iS 



Z 
Z 











' z 


z ' 




Q = 


i 


I 



s 



(6) 



Z 

z 



Then the matrix 



X k = (T- k T k ---T 2 T 1 )(°®W, 



is deterministic matrix factors away from the transfer matrix evolution T k - ■ -T\. However, 
it has the advantage that, unlike the product T k - ■ -Ti, it changes slowly as k varies. This 
can be seen from its evolution. From the definition we have 



(7) 



Note that in our case T k will be a small perturbation of T*. 

We call (Xk, < k < n) the regularized transfer matrix evolution, regularized at 
A*. In the next section, we will show that this evolution has a stochastic differential equation 
limit. But first let us check how to read off the eigenvalues of the operator M from X n . By 
(E]) the eigenvalues of M are given by 

A : det (((ri ® /2 ) Q ) n X n ) Q "(°"® 72 )) {lil) ) = 0} 

A : det (((Tj° ® J2)Q ) ri X n )^ 1 ) (1)1) ) = 
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Note that (T} 



{0®h)Q\n 



X„ is a matrix of the form 



Zn V 1711 y 

Zj JV\2 

JV 2 \ Zj v\ 22 

and so the (1, 1) block of Q(Ti°® l2)Q ) n X n Q- 1 is given by 

iS 11 (Z n+1 X 11 - Z n+l X l2 + Z n+1 X 2l - Z n+l X 22 )S ll = iSn{Z n+ \ Z n+1 )X n {I, -I) T S n 

Since Sn is nonsingular, have 

spec(M) = zeros A [det((Z n+1 , Z n+1 )X n (I , -I) T )] 

= zeros A [det %Z n+1 ((X n ) n - PQ12))] ■ (8) 

the second equality follows from the fact that the determinant is zero only for real A (as 
these are eigenvalues of a symmetric matrix M), and for real A, we have X 22 = X\\ and 
X 2 \ = Xi 2 . This will be shown in the proof of Proposition H] in the next section. 



4 The limiting of the transfer matrix evolution 



The goal of this section is to show that the regularized transfer matrix evolution introduced in 
the previous section has a stochastic differential equation limit. It will follow that the eigen- 
values of the corresponding operator converge in distribution to an explicitly constructible 
limit. 

We first state the scaling limit of the evolution of transfer matrices. 

Proposition 4 (Limiting evolution of transfer matrices). Fix r,a > 0. Consider the operator 
M n = rG x Z n + n~ l / 2 aV, where V = is a diagonal matrix with independent random 
entries of mean 0, variance 1 and uniformly bounded third moments. Denote the eigenval- 
ues of G by d±, . . . , d m and assume that < |A* — dj\ < 2 for all j and that the angles 
arccos((A* — rdj)/2), j = l,...,m are chaotic. Consider the regularized transfer matrix 
evolution corresponding to A*: (X^, < k <n). 

Then for any finite A C C we have convergence in distribution as n 00: 

(A$j(A* + X/n),t E [0, 1], A G A) (Y t (X),t e [0, 1], A e A) 

where Y is the solution of the SDE 



dY = S 2 



I 
-/ 



iXdtY + iaS 



dA 
-dB 



dB 
-dA 



SY, 



Yn 



(9) 
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and S is defined in (OJ). Moreover, for any fixed t, we have convergence in distribution of 
-> Y t (-) with respect to the uniform- compact topology of functions. The covariance 



X 



(n) , 
Lntj 1 



structure of the matrix-valued Brownian motions A(t), B(t) is as follows. We have 



A 



.1 1 ' 



Bjj = B 



(B is complex symmetric, not Hermitian) and 



E\Bt 



EAaAjj = EAaAjj E A 



jj 



Hj | 



(10) 



(ii) 



and all covariances that do not follow from the above are zero. 



Remark 5. In the case G = Z m by ([2]) we have 

1 



A(l) 



■■(A' + CI), 



B(l) 



1 



-.B' 



'12] 



y/m + 1 \Jm + 1 

Here A' is a version of the GUE with A\ ^ = Ajj standard complex normals for i ^ j and 
mean zero variance 1/2 i.i.d real normals in the diagonal. (The usual GUE would have 
variance 1 in the diagonal.) The random variable ( is standard normal and independent of 
A'. 

B' is an independent symmetric matrix with i.i.d. standard complex normal entries above 
and below the diagonal and i.i.d. mean zero complex normals with variance 3/2 in the 
diagonal (which are independent of everybody else). 

Next, we consider the eigenvalues of the operators M n . 

Corollary 6. Consider the operators M n of Proposition^ Assume that along some subse- 
quence Z n+1 — > Z as n — > oo. Then on this subsequence we have 



spec(M n ) zerosA det$S(Z(Y(l] 



Y(r 



(1,2) 



)) 



We proceed with the proof of Proposition HI The proof relies on the noise explosion 
phenomenon introduced in Section 1. 

Proof of Proposition^ We will assume a — 1, the general case may be proved the same 
way. We first study the evolution without scaling time. From (j7J 

X k - = (Tf^r^T- 1 - /)(WQ(Tj°^)* X k ^ 

The coefficient of above is given by 



Z k 

o z- k 



Q" 1 



o ±i + 4= vP 





Q 



z- k 
z k 



(13) 
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Note that the vector (X k (X) ■ X G A) is a Markov chain in k. In order to prove that it 
converges to the appropriate SDE we use Proposition [T3] of Section EJ 

We first rewrite Rk = -^R'k + w here R' k contains V k ° and R k contains XI from the 
middle term of Rk- We first focus on the noise term R' k which expands to 



" Z k 




I -z ' 




' 


v k °~ 




" z 


z ' 




' Z~ k 


o z- k 




-I z 












I 


I 




z k 



R' k = iS 

the three middle factors simplify and we get 



s 



iS 



z k v k °z- k 
z- k v k °z- k 



z k v k °z k 
-z- k v°z k 



s 



(14) 



(the diagonal block entries are negative conjugates of each other (note that V k is real) and so 
are the off-diagonals). For the convergence to the limit, we need to understand the covariance 
of the partial sums of Rk over k. (This is needed for condition ([25]) in Proposition [T3l .) 
For this, we may ignore the S factors for the moment, and study 



A i = Y^ Z k V k °Z~ k , B e = Y^ Z k V k °Z k 



(15) 



k=l 



The i, j entry of the first term is 



y^X Z i Z j) k ^ Olr V k,rOr,i = ^ ^ V k ,r{ Z i Z r ,jO r ,i 
k=l r r k=l 

Consider the complex covariance of the i,j and entries in Ag. Since that matrix is 
Hermitian, we may assume i < j and %' < j'. Since the v's are independent, the covariance 
is given by 



EAjjAj'j/ — ^^(ZjZjZj'Zj') k O r ,jOr,iOr,j'Or,i' ~ < 



k=l 



£(10^,10^) for 
im\\OA 2 ) for(z,0 = (j,iO 
(9(1) otherwise . 



This is because we can write ZiZjZi'Zy = e ia with \a\ < it and unless a = the sum will 
be 0(|a| _1 ). By the chaoticity condition a = can happen only if the conjugated are 
matched up in pairs to equal non-conjugated z i: in which case the indices have to be the 
same. The only options are i = = j' or i = j, i' = j'. 

This explains that the correlation between an entry A and the entry of B is always 0(1): 
there will be an odd number of conjugated z's so such a matching will not occur. 



10 



The expectation of the product on the other hand is given by 

i(OlOf,) for i = j and i' = f 



EAjjAj/ ji — ^ ^ (zjZjZjiZji) k 'y ^ O r jO rt iO r jiQ r ^ 



k=i 



0(1) 



otherwise . 



indeed, these are the only possible matchings for the conjugated and not-conjugated z's if 
i < j and i' < j'. Now 

EBijBi>ji = 0(1) 

since no matching can occur in the relevant product ZiZjZi/Zjr. Finally, for i < j and i' < j' 
(since B is a complex symmetric matrix) 



£(Ol Of) for i = %' and j = f 
0(1) otherwise . 



Turning to the drift term R'£ we first simplify it to get 

R>> = i(\ - \,)S 

This leads to the estimate 



/ Z 2k 
Z~ 2k -I 



t(A- A*) J2 S 



k=l 



I z 2k 
-z~ 2k -/ 



S = i(\ — \*)S 



s. 



tl 0(1) 

0(1) -u 



(16) 



s. 



;i7) 



Remark 7. From ( TH1) and ( TTBl it is clear that Rk is of the form 



a b 
b a 



and thus this will 



hold for / + Rk as well. The product of such matrices will have the same structure, which 
explains the last assertion of Section El 

We now turn back to the proof of Proposition HI For the proof of the convergence, we 
will use Proposition [131 f° r which we need to verify the conditions listed there. 



The boundedness of the cubic terms ff2^|) is proved as follows. It is clear that the 
increments X^+i — are bounded linear functions of A and the v i with coefficients 
given by the entries of X^. So condition holds as long as X^ is bounded. To 
ensure that Xk is bounded, we first consider a truncated process in which Xk will stop 
changing once \\Xk\\ > c. It follows from the proposition that the truncated process 
converges to the truncated version of the limit, for every c > 0. However, as c — > oo, 
the solution of the truncated SDE is with high probability equal to the non-truncated 
one. It follows that the processes converge. 
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Conditions ( )23|) and (I25|) in the truncated process follow from the calculations in the 
first part of the proof. Note a n (t,x) is a covariance matrix of dimension m 2 |A|. We 
demonstrate the process of checking its entries. First note that a n (t,x) is given by 



a n (k/n, x) 



E 



(^fe+i ~ x )ijA x k+i ~ x )i',f,y I x , 



(») 



X 



and the rows and columns are indexed by (i, j, A), (i', j', A'), respectively. Consider the 

A) = (i',j',X') case. 
Instead of working directly with a and b, we introduce 



n 



I Z 2k 
-Z~ 2k -I 



z k v£z- k 
-z- k v°z- k 



z k v£z k 
-z- k v°z k 



and let a n (t) and b n (t) be the n times the second and first moments of the m 2 |A|-vector 
y[ntj • Then we have 

Xk — x = iSY k Sx 

given Xk-\ = x. This is a linear function of x, and since x is bounded by truncation, 
it suffices to check fT2o^ for a and 6 rather than the original a, b. Similarly, ( l2l?|) for (the 
truncated) a, b is implied by 

|a n (t)| + |6 n (t)| < c 

This, in turn, is follows from the expression for Yj~ above (note that Z is diagonal with 
unit complex numbers). 

Returning to (125]) . we first check the existence of b so that 



sup 

t 



b n {s) ds 



b(s) ds 



0. 



This is clear by the following computation, based on (|T 



b(k/n) 



E4 W 



k \ 

^ n 



Z k IZ -k 

-z- k iz- k 



Z k IZ k 

_ z -k Iz k 



X- 

n 







+ 0(n- 1 ) 



I 

-/ 

where the O is uniform in k. To check the second moment terms, we first note that 
since the mean increments are of 0(l/n), it suffices to look at the variance d(t) instead 
of the second moment d(t). Then d(t) is the second moment of 



1 



z k v£z- 



Z kyO Z k 



-z- k v^z- k -z~ k v£z k 

This, by the computations above (starting at equation ( JT4l) ). converges to the covari- 
ance matrix of the corresponding noise term in the SDE. 
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For the last analyticity statement, we will need a bound of the form 

n\Xt\X)\\<f(X,k/n) (19) 

where f(X,t) is a locally bounded function. Then by Cauchy's theorem it follows that 
max^ ||x£ n) (A)|| is tight for any compact set A, and the claim then follows from Proposition 
[T4l We will follow the squared process instead. We have 

E(||X fc || 2 - ||X fc _i|| 2 | J- fc ) = Etr(AXXj_ l + AX T X k ^ + AXAX T ) 

n 

Taking expectations we now see that 



„ c(l + A 2 ) ., . 

Jl/ -Afc — Hi 1 1 yi. Jfe_l 1 1 ^ lJj ^Afe_i I 

n 



rearranging and taking a product for 1, . . . , k we immediately get the estimate 
E||X fc || 2 < (l + C(1 + A2) V llXoll 2 < exp(c(l + A 2 )t)||X || 2 



n 

from which the bound (1191) follows. □ 

Proof of Corollary By Skorokhod embedding, we can realize the distributional conver- 
gence of the matrix-valued analytic functions X^ =>■ Y(l) as almost sure convergence. In 
particular, along our subsequence, we have a.s. 

det((Z"+\ Z n+1 )X^(I, -if) det((Z*, Z*)Y X (1)(I, -if) 

uniformly on compacts, and the limit is analytic in A. Note that for A = the matrix in the 
determinant equals 

3(Z*(F (l)ii-n(l)i 2 ). 

whose distribution is absolutely continuous with respect to the distribution of the GOE; this 
follows from the SDE Thus the determinant as a function of A is not identically zero 
with probability 1, and the zeros of the subsequence converge to the zeros of the limit almost 
surely in our realization. The distributional convergence follows. □ 
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5 GOE as a limit 



In order to get the GOE limit, we need to take a sequence of boxes of the right dimension. 
First, some arithmetic conditions have to be satisfied, and by the following lemma they are 
in most cases. 

Lemma 8. Let A* be so that 

< 2 °> 

Then there exists Ti — > and rij — > oo so that for each n the points arccos((A* + ridj)/2) are 
chaotic and for each j , we have Hi arccos((A* + ridj)/2) — > mod 2tt. 

Proof. We will find rj so that the vector q = q{rj) = (arccos((A* + rjdj)/2) J=lr .. im has no 
nonzero integer vector orthogonal to it. It follows that the orbit {nq mod 2tt : n > 0} is 
dense on the m-torus, and we can find the desired convergence to 0. 

To find such rj we show that all but a countable set works. It suffices to show that for 
any fixed integer vector w, only countably many r« has q{ri) ■ w — 0. Note that r i— > q(r) ■ w 
is an analytic function on the interval [0,r*) where r* is the first value so that for some j 
we have Wj ^ and |A* + r*dj\ = 2. If this happens for a single j, then it follows that 
the function r i— > q(r) ■ w is not identically zero (since it has a singularity at r*), so it has 
countably many zeros. 

If this happens for more than one j, then the singularities may cancel out. However, this 
implies 

A* + r*dj = — A* + r*dji = 2, 
which is the reason we banished all such A*. □ 

Proposition 9. Let G be a real symmetric matrix with distinct eigenvalues di, and fix A* 
satisfying ( 120]) . Consider 

• a sequence of r v — > r so that the critical angles q v = (arccos(A* + r u dj))j=i y ,,^ m exist 
and are linearly independent over the integers, 

• a sequence n v so that for each j we have (n u + l)q u j — > mod 2ir and n v cha.(q u ) — > oo , 

• a sequence o v — > sufficiently slowly so that we still have maxj (e* (n ' y+1 - ) ' J '' J — l) — > 
and G y n u c\ia.(q u ) — > oo. 

• a sequence V v of diagonal perturbation matrices where the entries are independent, have 
mean 0, variance 1 and uniformly bounded third moment. 
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Then the regularized transfer matrix evolution for the operator 



n, 



rG x Z r 



v v - \J 



(21) 



satisfies the following. For any finite A C C we have convergence in distribution for the 
regularized transfer matrices 

{-{X[ n ^ (A* + \<T„/n v ) - I), t e [0, 1], A e A) (Y t (X), t e [0, 1], A E A) 

0~v 



where 



Y = iS 2 



' I 






A 


B 


Xt- 


-iS 






-/ 






-B 


-A 



s 



and S = limS^. Here A,B are matrix-valued Brownian motions with covariance structure 
given by (HDJ, (TTTj) . 

Moreover, the eigenvalue process of M v converges to the eigenvalues of 

S 11 Q(A(1) - B(1))S 11 . 

Remark 10. In the case when the graph is Z m , and S = si, the matrix 

S U $(A(1) - B(1))S U 

can be written as (m + l)~ 1 / 2 s 2 (K + bl), where b is a standard normal random variable, 
and if is a version of the GOE: a real symmetric matrix with independent mean zero real 
normal entries so that 

5/4 i = j 

Note that for the GOE, the diagonal terms have variance 2. The distribution of the matrix 
K is not invariant under orthogonal conjugation. The bl term amounts to a random shift of 
the local eigenvalue process. 

Proof. The proof follows very closely the proof of Proposition HI The only difference is 
that the frequency of oscillation in the oscillatory sums (fIBT) changes with v. In particular, 
for uniform convergence, we need to choose n u so that the geometric sums appearing in 
(fT5"|) are o(a u n u ). But these sums are bounded by cj cha,(q v ), which is still o{a v n u ) by our 
assumptions. 

For the convergence of eigenvalues, note that 



(22) 



L(z n+1 ,z n +- 1 )xto(i,-i) T 
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we write the above as 

-(Z n+l - 7, Z n+1 - t)X£\l, -if + - (J, 7).(/, -if + -1(7, /)(!? - /)(/, -If 

The first term converges to since X^\l, —If is tight and (Z n+1 — I) jo v — > by assump- 
tion. The second term vanishes, so considering the third term we conclude 

-(Z n+1 , Z n+1 )X^(I, -If (/, 7)F A (1)(J, -If 

all the convergence of the above terms are uniform of compacts. We now take determinants 
and use the fact that convergence of analytic functions implies the convergence of their zeros, 
the claim follows. □ 

Corollary 11. For almost all A* G [—2,2] there exists a sequence of integer 2-vectors 
(m,n) = (m u ,n u ) and r v — > 0,(7^,7^ so that the eigenvalue process of 

7 [(rZ m ) xZ n + aV- XJ] 

as v — oo converges locally to the Sine\ process. Here the noise matrices V u are diagonal 
with independent entries having arbitrary fixed mean distribution with uniformly bounded 
third moment. 



Proof. As n — > oo, by the methods of lErdos et al.l ( 12009bl ) the bulk eigenvalue process of the 



matrices K + bl (see Remark [TU!) converges locally to the Sinei process. The proof of this 
statement will be expanded in Lemma [16] in the Appendix. 
We now consider a A* so that 

2 — A* _^ cos(gi7r) 



2 + A* cos(g 2 7r) 

for any rational qi,q2- Consider a base graph G = Z m . Then by Proposition [9] and Lemma 
[8] we can find sequences n u , r v — > 0, a v so that the eigenvalues of ( 12T|) M n converge to those 
of s 2 (K + bl). After we have done this for a sequence of base graphs with m — > oo, we can 
find a diagonal sequence along which, after the appropriate scaling, the desired convergence 
holds. □ 

Remark 12 (Higher dimensional boxes). When G = Z mi x . . . x r L rrLd then the eigenvectors 
of G are of the form 

• • • ' kd) = n y^rr sin (^fr) 
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and the corresponding eigenvalues are 2 ^2 i=1 cos y t 

U(m + i)j (|OilMOil a > = II(i + ii 



I : Jf7r i . One can also calculate that 



kI=1 



If the eigenvalues are distinct and A* satisfies f[2"Ul) then Proposition will still apply. How- 
ever, the limiting symmetric real Gaussian matrix S (A(l) — B(l)) will have a more compli- 
cated covariance structure for d > 1 then the one described in Remark [10] and the current 
results are not strong enough to imply that the bulk scaling limit will be the Sinei process. 
One can check that the random matrix still contains an or der 1 component of a GOE 



Erdos et al. 



will go 



matrix which means that the local relaxation flow argument of 
through. The problem is caused by the fact that for the 'other' component the strong local 
semicircle result is not available. 



6 Appendix 

6.1 SDE limit of Markov chains 

The following is the main tool for pro ving convergence in the osciallat ory setting. It is 



Propo sition 23 in lValko and Viragi (120091 ). and is based on Theorem 7.4.1 of 
19861 ). 



Ethier and Kurtz 



Proposition 13. Fix T > 0, and for each n > 1 consider a Markov chain 

{X? G R d , £= 1... [nT\). 

Let Y™(x) be distributed as the increment X™ +1 — x given X" = x. We define 

b n (t, x) = nE[Y^ (x)}, a n (t, x) = nE[Y^ (x)Y^ (x) T ]. 

Suppose that as n — >■ oo we have 

\a n (t,x)-a n (t,y)\ + \b n (t,x)-b n (t,y)\ < c\x-y\ + o(l) 

sup£[|17(x)| 3 ] < en" 3 / 2 , 



(23) 
(24) 



and that there are functions a, b from R x [0, T] to M. d2 , M. d respectively with bounded first and 
second derivatives so that 

rt rt r t rt 

sup / a n (s,x)ds 

x,t Jq 









/ a(s,x)ds 


+ sup 


/ b n (s,x)ds 


lo 


x,t 


Jo 
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Assume also that the initial conditions converge weakly: 

Xq =>■ Xq. 

Then {XV" nt \ , < t < T) converges in law to the unique solution of the SDE 

dX = bdt + gdB, X(0)=X 

where B is a d- dimensional standard Brownian motion and g is any C 2 function satisfying 
gg T = a. 

We need a result that strengthens the convergence of random analytic functions. 

Proposition 14. Let f n be a sequence of random analytic functions on a domain D so that 
maxyi f n is tight for every compact A C D and f n =>■ / in the sense of finite dimensional 
distributions. Then f has a unique analytic version and f n =>- / in distribution with respect 
to local-uniform convergence. 

Proof. Pick a countable dense set of points D' C D. We can first find a subsequence so 
that the joint distribution of max^ \f n \ and f n (z),z G D' converges. By the Skorokhod's 
embedding theorem we can realize the sequence (f n (z), z G D') on a single probability space 
fl so that almost surely f n (z) — > f(z) for all z G D' and also max^ \f n \ converges (and so 
it is bounded) for A C A. By continuity, we can define the f n as analytic functions on D 
for the probability space Q. Then the f n have at least one analytic limit in the sense of 
uniform-on-compacts convergence by Montel's theorem. This limit must agree with / on all 
points z G D', so it is unique. 

By the above argument any sequence has a further subsequence that converges locally 
uniformly to some analytic / in distribution. But the distribution of / is determined by its 
finite dimensional distributions, so the limit is unique. □ 

Corollary 15. Let f n be random analytic functions on D so that ~Eh(\f n (z)\) < g(z) for 
some and increasing unbounded function h and locally bounded function g. Assume that 
fn{z) =>- f(z) in the sense of finite dimensional distributions. Then f has a unique analytic 
version and f n =^> / in distribution with respect to local-uniform convergence. 

Proof. By Cauchy's theorem maxA \f n \ is tight for every compact set A, and the claim follows 
from the Proposition. □ 
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6.2 The point process limit of the modified GOE 



Consider the n x n symmetric random matrix H n = n -1 / 2 ^ + bl) where b is a standard 
normal random variable, and K is the version of the GOE considered in ( 122]) . We will 
show that the eigenvalue process of H n in the bulk converges locally to the Sinei process, 
the bulk scaling limit of the GOE eigenvalue process. Bulk scaling means that we consider 



p(X)y/n(H n — XI) where |A| < 2 and p(X) = ^:\/4 — A 2 1|a|<2 is the semicircle density. 



Lemma 16. For any \X\ < 2 and compactly supported continuous test function O : 



we have 
da i 



da n Q(a u . . . , a k )p(X) "(p^ - PgoeJ [ x + 



,(*) 



Oil 



...A + 



Oik 



-►0. (26) 



np(Xy p(^) n , 

Here p$ n is the k-point intensity function of the eigenvalues of H n . The function Pgoeu 
is the same for the n x n GOE matrix with variances 2/n and 1/n on the diagonal and 
elsewhere, respectively. 



Proof. Since the proof works exactly the same way for any | A | < 2, we wil 



Our proof relies on the local re 
strong local semicircle law proved in 



axation flow argu ments of 



Erdos et al. 



Erdos et al. 



assume A = 0. 
f j2009bh and the 



( 120101 ). The argument in a nutshell is the fol- 



lowing: if we have a real symmetric random matrix whose eigenvalues are well-approximated 
by the semicircle law locally then by adding a small constant times an independent GOE 
matrix (i.e. by running Dyson's Brownian motion for a small time with our matrix as the 
initial condition) the scaled eigenvalue process of the resulting matrix will be close to the 
Sinei process. 

We will use Theorem 2.3 of lErdos et al.l (120101 ) (see also the comments after the theorem) 
which provides a powerful quantitative version of the previous argument. Since we can write 
-j^K — Ki + Ki where K\ is 1/2 times a GOE and K2 is an independent symmetric random 
matrix satisfying the conditions of the theorem then for any q > we have 

'75 dp,y/n 



2q 

x (/ } 



dai . . .da n Q(oti, . . . ,otk)p(0) h 



(27) 



K,n 



(fc) 

Pgoe,t 



" + ^)l^"- 1/1+e <«- 



-1/2} 



np(oy 

(We apply the theorem with E — 0, b — q/y/n,t = l/2,e' < 1 and 6=1— e' .) Since the 
spectrum of H n can be obtained by a random shift of the spectrum of n~ l / 2 K we have 



00 



npifS) ' ' np(0) 



db 



V2 



Tin 



-h±n (k) 

-- e 2 PKn 



Oil 



np(0) ' • • • ' + np (o) J ■ 



(2? 
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By Fubini for any nonnegative function F we have 



Using this with ( 1271) and (l2"8l) we get the upper bound 

/oo 2 
e-^\j-q 2 Cn- 1/i+£ (q- 1 + q~ l/2 ) < C'n~ l/A+£ 

which shows that the limit of the eigenvalue process of H n is the same as the limit of GOE 
eigenvalues multiplied by n and shifted by an independent Gaussian b n of variance n. We 
first condition on the sequence of these independent Gaussians, so the GOE eigenvalues are 
now centered at — b n , which, on the semicircle scale, converges to 0. In the classical literature, 
convergence is usually proved around a fixed w indow at c, |c| < 2 (on the semicircle scale). 



The moving window case is rigorously proved in IValko and Viragi (120091 ) . The lemma follows 



after removing the conditioning. □ 
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